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TIME  DEPENDENT  SOLIDIFICATION  OF  BINARY  MIXTURES* 


Bruno  A.  Boley 

Technological  Institute,  Northwestern  University 
Evanston,  Illinois 


Abstract 


The  problem  of  change  of  phase  in  a binary  mixture  is  considered  under  arbi- 
trary cooling  or  heating  on  the  surface.  Significantly  different  behavior  is 
noted  between  the  known  solution  for  a sudden  Jump  in  surface  temperature  and 
other  cooling  histories.  The  solution  is  presented  in  series  form,  and  a 
numerical  example  is  given  • *. 


Introduction 


The  problem  of  solidification  (or  melting)  of  mixtures  has  received  attention 
for  a number  of  years  both  from  fundamental  and  practical  points  of  view 
(cf.[l]).  Analytical  solutions  of  the  corresponding  coupled  heat  and  mass 
transfer  boundary-value  problem  have  been  discussed  for  semirinf inite  bodies 
under  sudden  changes  of  surface  temperature  (e .g. [2, 3 .Aj) . The  response  to 
these  special  conditions  has  been  found  to  be  characterized  by  similarity,  and 
by  constancy  of  concentration  in  the  solid,  and  of  both  phases  at  the  inter- 
face. This  behavior  cannot  be  expected  to  prevail  in  any  more  complex  pro- 
blem, whether  the  additional  complexity  be  due  to  different  cooling  histories, 
different  geometrical  configurations  or  other  effects.  The  purpose  of  the  pre- 
sent work  is  to  examine  the  problem  of  a semi-infinite  slab,  solidifying  under 
arbitrarily  time-dependent  cooling  conditions  at  the  surface,  and  to  derive  a 
solution  valid  for  short  times  after  the  start  of  solidification.  In  partic- 
ular, it  is  noted  that  the  special  behavior  alluded  to  earlier  is  not  exhib- 
ited in  the  present  solution. 


*Thls  work  was  supported  by  the  Office  of  Naval  Research 


Formulation  of  the  Problem 


Consider  a slab  (0<3t<L) , Initially  (t“0)  liquid  with  temperature  T and  with 

uO 

solute  concentration  C,  , under  presecrlbed  cooling  conditions  at  x>>0  and  (to 

Lo 

fix  ideas)  insulated  and  with  no  mass  transfer  at  x»L.  Solidification  will 

start  at  x-0  at  a time  t ,i.e.,  when  the  temperature  reaches  the  value 

m 


in  m 


T - mC, 
o Lo 


(1) 


mC  is  the  equation  of  the  liquidus  curve  in  the  (linearized) 

irf 


where  [2]  T » T 

o 

phase  diagram  of  the  binary  mixture  in  question.  The  pre-solidification  solu- 
tion is  easily  obtained  by  standard  methods  [Sj,  since  it  does  not  involve  any 
coupling  between  heat  and  mass  flow,  and  in  fact  corresponds  to 


C^(x,t) 


Lo 


tst 


(2) 


For  t>t  , the  solid  phase  will  occupy  the  space  0<5c<8(t),  and  the  problem  is 
m 

then  described  by  the  following  equations  (for  properties  constant  but  not 
necessarily  equal  in  the  solid  and  the  liquid): 

Field  equations: 

- 0 ; DCj^"  - - 0 in  s(t)  Oc<L 

Vs"  ■ ° " Cg  - 0 in  0<x<s(t) 

Boundary  conditions  (for  example): 

0 


-KgTg'(0,t)  - Q^(t)  or  Tg(0,t)  - Tg^(t);  Cg'(0,t) 

-K^T^'(L,t)  - Qj^(t)  or  Tj^(L,t)  - 


Interface  conditions: 


Tj^(s,t) 


Ts(8.t) 


T (t)  sT  - mC  (s.t);  C (s,t) 
m o L a 


kCj^(s,t) 


Vs'  ■ Vl’  ■ ’ ^^sS'  ‘ Vl'  “ V 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 


To  these,  initial  (t-t  ) conditions  must  be  added,  stipulating  that  temperature 

in 

and  concentration  are  continuous  at  t>t  and  that 

m 

8(t  ) - 0 (9) 

n 

In  the  above  equation  k is  the  partition  coefficient,  and  the  other  symbols 
have  obvious  meanings  (cf.[2]) 

Basic  Considerations 

It  is  convenient,  in  order  to  construct  a solution  to  eqs.  (3)  to  (9),  to  begin 
with  some  general  considerations  based  on  experience  with  the  analogous  un- 
coupled heat  conduction  problem.  This  eliminates  the  necessity  of  employing  a 
direct  awthod  of  solution  (e.g.,  the  embedding  technique[6])  %>hlch,  although 
applicable  to  the  coupled  problem,  is  likely  to  be  rather  cumbersome,  What 
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will  then  be  used  will  be  an  inverse  method  , that  is  one  in  which  the  form  of 
the  solution  is  assumed  at  the  outset,  and  it  Is  then  shown  that  all  conditions 
of  the  problem  can  be  satisfied  on  the  basis.  To  complete  the  solution  It  Is 
then  necessary  to  show  that  the  solution  thus  found  Is  the  only  possible  one. 
The  appropriate  uniqueness  theorem  will  be  presented  In  a subsequent  publi- 
cation and  is  not  considered  here. 


We  first  note  that,  in  the  uncoupled  problem,  a distinction  must  be  made  [7] 
depending  on  whether 

11m  %{y)ljy  - 0 or  lim  §(y)/v6^  i 0 


y>0 


y-»0 


where 


y - (t/t  ) 
in 


(lOa) 


The  second  of  (10)  is  satisfied  by  the  similarity  solution  earlier  mentioned, 
which  is  a direct  extension  of  the  classical  Neumann  solution  [5]  of  Stefan's 
problem.  We  note  that,  in  both  uncoupled  and  coupled  solutions  of  this  type, 
the  form  of  § is  the  same,  i.e.,  § = for  y^l,  where  X is  a constant.  We 

will  return  to  that  problem  later;  our  present  aim  is  the  discussion  of  prob- 
lems characterized  by  the  first  of  (10),  and  we  shall  assume  that  there  too 
the  form  of  ?(y)  is  unchanged  by  coupling. 


The  form  of  the  function  ?(y)  is  known  in  a number  of  cases.  The  starting 
solution  (i.e.,  the  first  term  of  series  expansion)  is  know  for  any  arbitrary 

cooling  history  [7].  For  the  general  class  of  problem  in  which  the  surface 
heat  flux  is  expressible  in  a series  of  a half-integral  powers  of  y,  it  is 
found  [7]  that  5(y)  is  also  so  expressible,  while  for  the  companion  problem  in 
which  the  surface  temperature  is  so  expressible,  5(y)  requires  [8]  a series  in 
powers  of  y^.  Let  us  consider  the  former  of  these  classes  of  problems,  and 
assume  that  the  heat  flux  is  continuous;  then  [7]  gives 

?(y)  - \ (11) 

n"0 

Inspection  of  the  second  of  eq.(8)  can  now  be  used  to  guess  the  first  term 
of  a series  representation  for  the  concentration.  The  right-hand  side  of 
that  equation  is  proportional  to  s,  or,  with  (11),  initially  proportional  to 
Hence  It  Is  reasonable  to  expect  that  the  concentration  be  distributed, 
for  short  times.  In  the  isanncr  corresponding  to  a surface  flux  also  proportional 
to  or 


3- 


(12) 


C_  (X ,y ) 2 

— — = 1 + £^2!^^  erfc(XB^6/v^)  X^l 

bO 


where 


" * ^ * \.s  “ Ks\,s  • ^ 

a m 


(13) 


For  C we  write  a similar  expression,  but  imaged  about  X • 0 so  as  to  satisfy 
S 

(5),  or,  still  for  y<l. 


C_(X,y)  f 7 2 “I 

— ? * + — — yj^i  erfc(XBg/^/91  + i erfc (-XBg/v^ J 

'"Lo 


(14) 


where  the  last  of  (7)  requires  that  Cg^  “ *^^Lo’  where  the  factor  k/2  has 
been  Introduced  for  future  convenience. 


To  extend  the  solution,  it  is  again  conjectured  that  a series  in  half-integral 

powers  of  y is  appropriate,  or 

00 


C-(X,y)  ^ n/2n 

^ + Z.  ^^Ln^  ^ erfc(XBJ^6/^/9■  ) 

Lo 

C-(X,y) 


n=2 


"Lo 


k|l  + Cgjjy"'^^[i”erfc(xBg/^^)fl‘'erfc(^-XBg/v9)1} 


(15a) 


(15b) 


n=2 


Similar  argruments  may  now  be  applied  to  the  first  of  eqs.(8)  to  obtain  the 

form  of  the  temperature  functions.  In  this  case,  one  must  be  careful  however 

to  add  particular  solutions  of  (3)  and  (4)  which  will  insure  satisfaction  of 

the  non-homogeneous  boundary  conditions  (5).  The  simplest  way  of  achieving 

. * 

this  is  to  employ  [6,9]  the  analytic  continuation  T (X,y),  into  the  post- 
solidification period,  of  the  pre-solidification  solution.  Thus 


Vj^(X,y)  = 


Vg(X,y) 


:L(x.y) 

* " 

T (X.y)  , y, 

T-  -T  (t  ) L 
Lo  n IQ  A 

n-2 

TjW.y) 

_ T*(x.y)  , 

(t  ) 

(t  ) 

(16a) 


^E^Tsny"^^[i"erfc(X/^  + i"erfc(-X/v9)] 


n«2 


(16b) 


Eqs.(15)  and  (16)  satisfy  all  conditions  of  the  problem,  with  the  exception  of 
the  conditions  to  be  satisfied  at  the  solid-liquid  Interface,  namely  (7)  and 
(8).  In  dlswnslonless  form,  these  take  the  following  form  (the  arguments 
^,y  being  understood  throughout): 


4 


(17a, b) 


c , 

< 

II 

< 

= V ^1  - ' 

L S 

^Lo^ 

(C„/C 

= k(C  /C 

\ S Lo/ 

\ L Lo/ 

ax 

ax  M dy 

) d C /C.  I 

S Lo 

/ f S)  \ L Lo/ 

ax 

\h  ' ax 

IS 

dy 


(17c) 

(17d) 

(17e) 


where  the  dimensionless  notation: 

2 I ’ T 


mC. 


M 


Lo 


' "o=  T 


In  ’ 

LO  m 


(18) 


!l 

" “ ^ ’ A ' ' 

The  five  sets  of  coefficients  S_,C,  ,C  ,T  ,T  must  be  adjusted  so  as  to 

Ti  Ln  Sn  Ln  Sn 

satisfy  the  five  equations  (17a-e). 


Solution  and  Discussion 


The  solution  now  requires  the  substitution  of  (11),  (15)  and  (16)  into 
(17a-e);  separation  of  the  resulting  equations  into  terms  of  like  powers  of  y, 
then  gives  sets  of  equations  from  which  the  unknown  coefficients  can  be  derived. 
The  process  is  straight  forward  but  rather  laborious,  since  it  requires  expan- 
sion in  powers  of  y of  the  several  Integrated  error  functions,  whose  arguments 
themselves  are  power  series.  The  details  of  the  process  will  be  omitted,  but 
it  may  be  noted  that  the  calculations  fall  in  a distinct  pattern.  It  is  con- 
venient to  consider  first  the  highest-order  terms  in  eqs.  (17a, d,e);  these 
yield  three  simultaneous  equations  for  (?o*^L2*^L2^  * coefficients 

^^S2'^S2^  are  then  found  from  (17b, c)  respectively.  The  pattern  repeats  for 
further  coefficients:  a set  of  three  simultaneous  equations  for  (^j^, 
results  from  the  next-order  terms  in  (17a, d,e),  with  Cl^s3»^g3)  obtained  again 
from  (17b. c),  and  so  forth. 


As  an  example,  the  case  of  a constant  flux  Q at  x>0  was  considered  in  detail. 
Here 


I 

T 


*(x.y) 

Lo”'m<^m> 


Lo 


T -T  (t  ) 
o m m 


- yrrfl+y)  ierfc 


7^ 


(19) 


5 


and  the  following  results  are  obtained  (for  simplicity  we  take  6 * T)  * 1 and 


\ ■ ‘s>‘ 

First-order  terms; 


? = i 

° Is  + 3^  B(l-k)AV^ 


C,,  = C,,  = 6^  B(l-k)?^  ; T - = T , = 3tt?  /M 


(20a) 


Second-order  terms; 


5 

^ + 16(l-k)BAV 

M O 


■ S3  ■ 32(1-M8!j  ; T^3  ■ T^3  - lVS-?i/« 

Third-order  terms: 


(20b) 


4a^g^^+2M(3^g^-l)+3TT^^^g^(^  -•3g>3BVM(l-t^)g![32-15TT(l-k)] 


+2(l-k)BAV^M  ’o'  o 


^ ' 15tt  + 307TrAV^BM(l-k) 

= 30^(l-k)B^252  + 3^B(l-k)§^^]  ; - 1928^ (1-k) 


(20c) 


^L4  ■ ^'52  ^ % - %%]  ■ hi  ■ V - 

and  so  forth.  Some  comments  on  the  above  solution  may  be  useful. 

The  method  employed  for  the  solution  makes  It  clear  why  a distinction  must  be 
made  depending  on  which  of  eqs.(lO)  holds.  In  fact,  when  (as  at  present)  the 
first  of  (10)  is  valid,  the  various  error-function  Integrals  tend,  as  y-»0  and 
on  X - §(y)*  to  constants  Independent  of  the  problem  parameters.  When  the 
second  of  (10)  holds.  In  contrast,  the  error  functions  are  dependent  on  the 
proportionality  constant  X between  5 ®nd  a/V • This  Implies,  for  example,  that 
Cj^(§»y)-»f  (X)  as  y-*0,  and  therefore  T^(t)  la  also  dependent  on  X.  The  quantity 
X would  appear  in  all  equations,  and  the  present  separation  of  the  several  equa- 
tions would  no  longer  be  possible.  The  procedure  followed  here  could  of  course 
still  be  employed,  but  It  would  be  much  more  cumbersome.  It  Is  rather  more  con- 
venient to  solve  the  problem  separately  In  these  cases,  as  Indeed  was  done  In 
[2.3,4]. 

The  present  solution  reduces  to  the  uncoupled  one  [9]  if  either  A«o  or  k-1;  In 
all  other  cases  the  form  of  ^ is  the  same  but  the  coefficients  differ.  The 
concentration  of  the  Interface,  given  by 


1 

k C, 


1 + ^ 


3/2 


n^-7^L2,'’y  * -(21) 


So  “ “Lo 

exhibits (Fig. 1) a gradual  Increase  from  Its  initial  value.  Note  that  the  concen- 


tration In  the  solid  Is  not  uniform. 
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LIST  OF  SYMBOLS 


Symbols  listed  in  the  [brackets]  are  dimensionless,  and  are  defined  in  the 
equations  indicated. 


[A,  eq.(18)] parameter  characterizing  the  liquidus  line 

[B,  eq.  (13)] * 

c specific  heat 

C mass-fraction  of  solute 

D diffusion  coefficient 

f function 

k partition  coefficient 

K thermal  conductivity 

I latent  heat  of  fusion 

L thickness 

m slope  of  liquidus  line 

[M,  eq.(18)] principal  fusion  parameter 

Q heat  flux 

s,  [^,  eq.(lOa)] position  of  solid-liquid  Interface 

[y<  CQ*  (10a)]... time 

T,  [V,  eq8(16)] temperature 

X,  [X,  eq.(13)] distance  from  exposed  surface 

[6,  eq.(13)] 

[ti,  eq.(18)] ratio  of  liquid  to  solid  dif fusivities 

K dif  fuslvity 

[X] similarity  coefficent  in  constant 

temperature  solution 
p density 


Fig.l.  Variation  of  interface  concentration  Cj^(?,y) 
interface  position  5(y)  with  time. 


0.1;  k ■ A “ 0.5;  T|  ■ 6 ■ B 


B. 


2Cg(5,y)  and  of 
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